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ABSTRACT 

An approximate equation is derived for calculation of gas weight 
flow rate through subsonic orifices for small differential pressures. 
Comparison to the exact equation shows good correlation even at con- 
ditions far removed from those assumed in the derivation. A limited 
test program verifies the applicability of the approximation. 
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SYMBOLS 


A 12 = Orifice cross section area 
P-^ = Pressure upstream of orifice 

P 2 = Pressure downstream of orifice (thruster inlet pressure) 

fcp - *1 - p 2 

= Temperature upstream of orifice 
g = Acceleration of gravity 
Cd - Coefficient of discharge 
R “ Gas constant 
CO - Exact flow rate offgas 
(jJ a = Approximation for VO given by eq. (7) 

W m = Measured flow rate of gas 
y = Specific heat ratio 


INTRODUCTION 

The use of sonic orifices in the measurement of gas flow rates 
between 10”5 and 10”4 lbs/sec for laboratory performance evaluation of 
millipound thrust ammonia resistance jet thrusters has been common 
practice (Ref.'s (1) and (2)). This type of testing, however, usually 
does not fully simulate the operation of the thruster on a spacecraft. 

The prime requirement on a spacecraft is for a constant thrust which 
dictates that the inlet pressure to the thruster be held constant. For 
a long gas pulse, the specific impulse (thrust/^) of an ammonia resistance 
jet could vary by a factor of two from initiation to completion of the 
pulse due to the decrease in gas temperature at the thruster nozzle as 
the temperature of the resistance jet drops. This means that if the thrust 
is held constant, the flow rate at the completion of the pulse would be 
double the initial flow rate. If a sonic measuring orifice is placed in 
front of the thruster, the requirement of a constant downstream pressure 
(P 2 ) is very difficult to maintain because Pj has to be continuously varied 
as the flow rate increases. With a subsonic orifice, however, P 2 can be 
held essentially constant with only small variations of P]_. 

APPROXIMATION OF ORIFICE FLOW EQUATION 


The exact expression for the weight flow of a compressible fluid 
through an orifice is given by 



eq.(l) 


By using a series expansion for the pressure ratio terms under the radical 
in eq. (1), Ref. (3) derives an approximation for the orifice flow equation 
when AP is small: 
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neglecting terms in (■ A -£ ) to the third power or greater, the following 


By 

equation is obtained: 
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eq.(4) 


Substitution of eq.(4) into eq.(l) yields 

Oj-V A .12 c D (1 . _3_ Ai?j- 

Since eq.(5) is valid only for small AP's, 
assume that the term (3 /2 )( AP/p^) should also be small. Compressible 
gases of interest have specific heat ratios of 1.0 to 1.67. It is reason- 
able to further assume, therefore, that setting (3/2)0 *#1 would introduce 
negligible error into the flow rate approximation. Eq.(5) thus becomes 


eq. (5) 

it is reasonable to 
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eq. (6) 


For the application of interest, the gas temperature at the orifice 
is essentially constant at ambient conditions. If one further assumes 
that Cp = constant for a particular orifice, eq.(6) can be written 


Kg (P 2 AP) 0 ’ 5 


eq. (7) 


where 05 

K g = a 12 c D < — * • e£ l-< 8 > 

Thus with Ap small, y ££1.5, Cp = constant and = cons tan t^ the weight 
flow varies linearly with (P 2 A p )^*^ e 


To test the validity of the assumptions used in deriving eq.(7) 
calculations were made on a digital computer comparing the weight flow 
calculated by eq.(7) to that calculated by eq.(l). Various gammas 
were used and the pressure ratio was varied from P 2 /Pl =1.0 down to the 
critical pressure ratio. Precentage error was found from 

% Error = IM z MJ olI x 100 

ijj . eq. (9) 

It can be shown that percentage error is independent of Ai 2 , R, Cd, g 
and Tit and can be expressed as the following function of (P 2 /Pi) and jj' : 
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x 100 




eq.(10) 


Table 1 summarizes the data found from this analysis. 

TABLE 1: % Error = x 100 

00 



1.14 

1.31 

1.4 

1.5 

1.66 

.98 

.32 

.15 

.07 

.0008 

.10 

.95 

.83 

.38 

.19 

.005 

.24 

.90 

1.72 

.79 

.40 

,02 

.49 

.85 

2.68 

1.24 

.63 

.005 

.73 

.80 

3.73 

1.73 

.90 

.09 

.98 

.75 

4.88 

2.27 

1.19 

.15 

1.23 

.70 

6.14 

2.88 

1.53 

.24 

1.48 

.65 

7.53 

3.55 

1.91 

.34 

1.74 

.60 

9.08 

4.31 

2.35 

.48 

1.99 

(P2/Pl) crit=.577 
(P 2 /Pl.) crit=.544 
(P 2 /P 1 ) crit=.528 
(P 2 /P 1 ) crit=.512 
(P2/Pl) crit=.488 

9.85 

5.28 

3.10 

.83 

2.53 


Additional calculations show that using eq.(5) introduces less 
error than the use of eq. (7). This is to be expected because eq.(5) 
uses the correct value of y . Eq.(7), however, expresses uj as a linear 

function of (P2AP)®»5, and is more convenient when calibrating an 
orifice. 

* 

The data of Table 1 shows that the approximate equation for 6C* gives 
good correlation with the exact equation even at conditions far removed 
from those assumed in deriving eq.(7). As expected the amount of error 
decreases as t \ P and approach zero and 1.5 respectively. Table 1 
shows, however, that even at critical flow, the error introduced by 
using eq.(7) is less than 10% for a gas with 1.37 and approximately 

2.5% for a gas with K = 1.66. Eq.(7) gives a high estimate of the flow 

rate for the four gammas considered with X £ 1.5 and a low estimate for 
the y = 1.66 gas . 


The principal value of eq.(7) and the limited supporting orifice 
calibration data, besides providing a simpler equation than that of eq.(l) 
for estimating the flow rate through an orifice, is the indication that 
only a few points need be taken in the calibration of the orifice to 
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obtain the straight line relationship of (jJ vs. (P 2 AP)^*"*. 

♦ 

An interesting observation is that a linear relationship between lO 
and (PiAP )* 5 can be obtained by assuming that (3/2 Y)/AP/P*j^| is 
negligible with respect to unity in eq.(5). This approximation, however, 
is less accurate (in comparison to the exact^ expression f or VO) than eq. 

(7) at all pressure ratios for gases with o^l.5. For gases with 
V > 1.5, this approximation gives more accurate results for very small 
AP's, but it gives significantly larger errors as the orifice pressure 
ratio approaches critical conditions (e.g. at y= 5/3, and P 2 /P 1 = .488, 
this approximation gives 40% error) . 

ORIFICE CALIBRATION 

To further verify that Vt) can be approxi.mated as a linear function 
of (P 2 AP)° °5, a number of small diameter orifices were calibrated with 
NH 3 . The calibration procedure is detailed in the Appendix. The 
results of the # calibration, shown in Figure 1, verify the linear relation- 
ship between 4 ^ and (P 2 AP)®*^* The values obtained for the discharge 
coefficients for these orifices, which were made by drilling holes in a 
sheet of .0625 inch thick stainless steel, varied from .77 to .94 (assum- 
ing the orifice diameter measurements were correct) but appeared to be 
constant for each particular orifice in the range tested. 

The orifice calibrations were limited to cases where the pressure 
ratios (P 2 /P 1 ) were large (the calibration was done to satisfy the flow 
requirements of a resistance jet test program) . More extensive testing 
would be required to verify that Cj), for a particular orifice, is constant 
over a larger pressure ratio variation. 
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APPENDIX 


ORIFICE CALIBRATION PROCEDURE 

Figure A shows the schematic of the orifice calibration system. The 
orifice is calibrated by the following procedure: 

(1) Fill the sphere with gas. 

(2) Weigh the filled sphere. (use 10,000 gram capacity analytic balance 
with a jp.l gram readability). 

(3) Attach the sphere to the flow system as shown in Figure A. The 
pressure regulator and valve (B) are preset by trial runs to yield 
the desired P 2 andAP. 

(4) Open valve (A) and start a stop watch at the same time. Let the 
gas flow for a sufficient amount of time (depending on orifice size) 
to offset transient and measurement errors. 

(5) Close valve (A) and turn off stop watch at same time. 

(6) Remove the sphere from the flow system and weigh it. 

(7) The flow rate is found by subtracting, the final weight of the 

sphere from the initial weight and dividing by the time period for 
which the gas was flowing. 9 

(8) Vary Pj, and AP over a large enough range to plot the curve of CO 
vs (P 2 AP)®*^ (should be a straight line for small AP's) for the 
range of interest. 

The calibration system shown in Figure A exhausts the gas leaving 
valve (B) at ambient pressure. This valve can be replaced by an orifice 
which exhausts to a vacuum (to simulate the thruster). Some orifices 
were calibrated by both techniques and gave identical calibration curves. 

The water bath shown in Figure A is used to maintain the test gas 
at a uniform temperature by means of a temperature controller. 




FIGURE A: ORIFICE CALIBRATION SYSTEM 
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